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We generalize the (2 + l)-dimensional Yang-Mills theory to an anisotropic form with 
two gauge coupling constants e and e'. In axial gauge, a regularized version of the 
Hamiltonian of this gauge theory is H + e' 2 Hi, where H is the Hamiltonian of a set 
of (1 + l)-dimensional principal chiral sigma nonlinear models and Hi couples charge 
densities of these sigma models. We treat Hi as the interaction Hamiltonian. For gauge 
group SU(2), we use form factors of the currents of the principal chiral sigma models 
to compute the string tension for small e', after reviewing exact S-matrix and form- 
factor methods. In the anisotropic regime, the dependence of the string tension on the 
coupling constants is not in accord with generally-accepted dimensional arguments. 
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1 Introduction 



In this paper we calculate the string tension of pure (2 + l)-dimensional Yang-Mills 
theory at weak coupling. This is done in an anisotropic version of the theory, with two 
different, but small, coupling constants. The result clearly establishes that confinement 
is not restricted to the strong-coupling region. 

The method used is a development of the proposal of reference PQ, which in turn 
was inspired by Mandelstam's arguments for confinement in QCD [2]. The starting 
point is a regularized Hamiltonian in an axial gauge, A± = 0, where Aj, j = 1, 2 is the 
Lie-algebra valued gauge field (the component Aq is also set to zero). The Hamiltonian 
may be written as a sum of two terms, namely 

and 



e 2 



H 1 = -j d'xTtSt, 

where Sj, j — 1, 2 are the components of the electric field and B is the magnetic field 
B = i[di— L4i, d 2 — L4 2 ] = diA 2 Though we may write S 2 = —iS/5A 2 , in the Schrodinger 
representation, the formula for Si is a nonlocal expression, obtained from Gauss's law 

m 



Si{x) = -/ dy l [d 2 -iA 2 {y\x 2 ),S 2 {y\x 2 )] 
i 

dy'V^^x^-S^y 1 ^ 2 ), (1.1) 



where T> 2 is the adjoint covariant derivative in the two-direction. The local form of 
Gauss's law is explicitly satisfied with (jl.lj) . provided a residual gauge invariance 

dx l V 2 S 2 ^ = , (1.2) 

is satisfied by physical states \l/ (this condition must be modified slightly if quarks are 
present). 

If the theory is regularized on a lattice, H is the Hamiltonian for a set of decoupled 
SU(N)xSU(N) principal chiral nonlinear sigma models. These sigma models are coupled 
together by the interaction Hamiltonian Hi. We shall treat the coefficient of Hi as 
small, meaning that we consider an anisotropic modification, 



Hl = ^J d 2 xTrS 2 . 
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By assuming e! -C e, simple arguments were made that the theory is in a confining 
phase [I]. It was also suggested that Rayleigh-Schrodinger perturbation theory in e' 
was infrared finite - but it is now clear to the author that this is not the case. A better 
understanding is needed to obtain quantitative results. 

In this paper, we shall work with gauge group SU(2) and exploit exact knowledge 
of the form factors of the (1 + l)-dimensional 0(4) ~ SU(2) xSU(2) nonlinear sigma 
model [3 . The expressions for the form factors are valid for small e, as we shall show, 
so both gauge couplings must be small. After reintroducing the component Aq of the 
gauge field, we integrate out A 2 using these form factors. We determine the resulting 
effective action of A to leading order. This effective action is infrared finite. 

We shall find a result for the string tension in the x 1 -direction which is different 
from that in reference PQ and will discuss how the difference comes about. We do not 
calculate the string tension in the x 2 -direction here, which is still under investigation. 

In the course of these investigations, we discovered a paper by Griffin which suggests 
a program very similar to ours [3]. Evidently, his ideas were taken no further. Griffin's 
proposal was for the light-cone gauge, instead of the axial gauge, but the two are very 
similar in certain respects. 

We should say a few words about other analytic approaches to Yang-Mills theories 
in (2 + l)-dimensions. Probably the best known paper is that of Feynman [B]. Though 
some of the mathematical details of Feynman's paper are not correct, evidence for some 
of his conjectures were found in reference j^J, which utilized some methods introduced 
in reference |7j. In particular, the low-magnetic-energy region of configuration space 
was argued to be bounded in (2 + 1) dimensions (it is unbounded in (3 + 1) dimensions 
jjj), which is indicative of a mass gap. 

An approach, which at first appears rather different, was developed by Karabali and 
Nair [S] in which new coordinates for configuration space were devised. A re-derivation 
of one form of the Hamiltonian discussed in reference [S] was done in [5| using a kind of 
non-Abelian exterior differentiation; in this way a connection was made with the ideas 
of references j7] and jHj. 

The major achievement of reference [S] was the finding of a mass gap at strong 
coupling and a resummation of the strong-coupling expansion. This resummation was 
used to obtain a string tension which is independent of the ultraviolet cut-off. We 
should mention, however, that very similar results can be obtained by analytic lat- 
tice methods. The strong-coupling spectrum in the lattice Hamiltonian formalism has 
the same dependence on the continuum coupling as in Karabali and Nair's formal- 
ism. A different resummation method on the lattice, due to Greensite JTU] also yields 
the strong-coupling vacuum wave functional. The form of Greensite's vacuum wave 
functional in the continuum limit in 2 + 1 dimensions is 

*-=^-i(irW/ toBi(ir 

The wave functional that Karabali and Nair obtained has the form, in the infrared 
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limit, 



C A e 4 



TV 



I 



d 2 xTr[Fij(x)] 2 



where Ca is the quadratic Casimir for the adjoint representation. Both vacua yield a 
string tension proportional to e 4 . It does not seem to have been noticed before that 
the essential features of Karabali and Nair's strong-coupling expansion and those of 
the lattice strong-coupling expansion worked out by Greensite are the same. 

To obtain a genuine proof of confinement though strong-coupling methods, it must 
be shown that the strong-coupling expansion has a finite radius of convergence in the 
1/e or an infinte radius of convergence in the dimensionless coupling l/po- This has 
not been accomplished yet. 

Recently a new approach has appeared [llj, which uses a scheme to improve a 
Gaussian Ansatz for the vacuum and excited-state wave functionals in the variables of 
reference jH]. 

What distinguishes the approach of pQ and this paper from other work is that does 
not exploit strong-coupling approximations or any Ansatz for wave functionals. We do 
make assumptions, which are quite different from the assumption of reference We 
assume that an anisotropic weak-coupling expansion makes sense and that the result 
for the SU(2)xSU(2) principal-chiral-sigma-model exact form factor [3| (which has been 
checked in the 1/N expansion for the 0(4) formulation) is correct. We think that our 
approach leaves no doubt that confinement occurs at weak coupling, granting that the 
anisotropy is something we would like to get beyond. Or perhaps not - the anisotropic 
theory is asymptotically free and not finite like the isotropic theory. In this respect it 
is more like real QCD in 3+1 dimensions. 

By treating the coefficient of Tr£ 2 (instead of TtB 2 ) as small, our method is in- 
herently a weak-coupling approach. We think much insight can be gained by working 
systematically at small coupling. Furthermore, to solve the much harder problem of 
QCD in (3 + 1) dimensions, a weak-coupling understanding is essential. 

There is a well-accepted argument concerning how the mass gap M and string 
tension a depend on the coupling constant (see for example references [H] and [T2]). 
We find that such an argument fails in the anisotropic weak-coupling regime. The 
argument goes as follows: Yang-Mills theory is a perturbatively-ultraviolet-finite field 
theory (there are severe infrared divergences, but let us ignore these). Thus, after being 
suitably regularized, the coupling has a nonzero finite value, as the ultraviolet regulator 
is removed. Since the coupling squared e 2 , has units of cm -1 , the mass gap must behave 
as M ~ e 2 and the string tension must behave as a ~ e 4 . In our anisotropic case, the 
same dimensional reasoning implies 



for some set of numbers p and P and dimensionless constants C p and Kp. Our final 
answer for the string tension is quite different from (jl.3j) . We show in Section 6 that 





p 



p 
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for two quarks separated in the x 1 -direction, 

e 2 47r 
cr ~ — exp — , 

a e A a 

where a is a short- distance cut-off. We believe that for e — e' the dimensional argument 
does yield the right answer, but that there is a crossover phenomenon between ()Oj) to 
the behavior we find in the anisotropic regime. 

Another application of exact form factors to the (2 + l)-dimensional SU(2) gauge 
theory has just appeared [T^]. In this work, the form factors of the two-dimensional 
Ising model are used to find the profile of the electric string, near the high-temperature 
deconfining transition, assuming the Svetitsky-Yaffe conjecture. 

An interesting question is the value of k-string tensions for gauge group SU(iV) (see 
reference [TJ] for detailed review of this matter). This is not discussed in this paper, 
since the gauge group is SU(2). Thus the value of k is always one. Another question 
is whether adjoint sources are confined. Both of these issues are addressed in a new 
paper [T^]. The sine law is clearly seen for the vertical string tensions. The situation 
is less clear for horizontal string tensions; at zeroth order in g' Q there is a Casimir law, 
but there are corrections. We are unable to calculate these corrections, because we do 
not know the form factors for SU(iV) x SU(iV) principal chiral sigma models. Adjoint 
sources are shown not to be confined. 

In the next section, we discuss the regularized Hamiltonian. In Section 3, we go to 
the axial gauge and split this Hamiltonian into the Hamiltonians of (1 + l)-dimensional 
0(4) ~ SU(2) xSU(2) nonlinear sigma models H and a nonlocal term Hi. We discuss 
how to find the effective action of the temporal gauge field in terms of correlators of 
the nonlinear sigma model in Section 4. In Section 5, we determine the leading-order 
effective action using the exact form factors of the 0(N) nonlinear sigma model in 
(1 + l)-dimensions. The static potential is then found between two quarks separated in 
the ar-direction in Section 6. The physical picture of confinement of glueball excitations 
and quarks separated in the x 2 -direction is presented in Section 7. We discuss some 
future endeavors in Section 8. We give a review for non-experts on the exact S-matrix 
[TB] and form factors j3] of the (1 + l)-dimensional 0(N) nonlinear sigma model in the 
Appendix. 

2 The regularized Hamiltonian 

We will quickly review the Kogut-Susskind Hamiltonian formulation of lattice gauge 
theory. If the reader finds this discussion incomplete, we refer him or her to the book 
by Creutz [T7j. 

Consider a lattice of sites x of size L l x L 2 , with sites x whose coordinates are x 1 
and x 2 . We require that x 1 /a and x 2 /a are integers, where a is the lattice spacing. 
There are 2 space directions, labeled j = 1,2. Each link is a pair x, j, and joins the 
site x to x + ja, where j is a unit vector in the j th direction. 
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We use generators tb, b = 1,2,3, of the Lie algebra of SU(2), which are related to 
the Pauli matrices by t b = a b /^/2. The identity matrix will be denoted by 11. 

For now, the Hamiltonian lattice gauge theory will be in the temporal gauge A = 0. 
The basic degrees of freedom, before any further gauge fixing, are elements of the group 
SU(2) in the fundamental (2 x 2)-dimensional matrix representation Uj(x) G SU(2) at 
each link x, j. The relation between these variables and the continuum gauge field 
is Uj(x) = e~ mAj ( x \ There are three self- adjoint electric-field operators at each link 
lj(x)b, b = 1, 2, 3. The commutation relations on the lattice are 

[lj(x)b,lk(y)c] = iV25 xy 5 jk e bcd lj(x) d , 



[lj(x) b , U k (y)} = -5 xy 5 jk t b Uj(x) , (2.1) 

all others zero. 

The Hamiltonian is 

^ = EEE |[^f-Ei [TrU 12 (x)+TrU 21 (x)} , (2.2) 

X j = l b=l x ^° 

where 

U jk (x) = Uj(x)U k (x + ja)Uj(x + ka) ] U k (x) ] , 

where j and k are the unit vectors in the j- and /c-directions, respectively, and the bare 
coupling constant go is dimensionless. The coefficient of the kinetic term is just half 
the square of the continuum coupling constant e, namely (?q/(2o) = e 2 /2. This is why 
the mass gap in (2 + 1) dimensions is proportional to e 2 in strong-coupling expansions. 

We denote the adjoint representation of the SU(2) gauge field by TZj(x). The precise 
definition is 1Zj(x) b c t c = Uj{x)tbUj{x)^ . The matrix IZj lies in the group SO(3). 

Color charge operators q(x)b, may be placed at lattice sites. These obey 

[q(x) b ,q(y) c ]=iV2e bca 5 xy q(x) a . (2.3) 
In the presence of static charges, Gauss's law is the condition on physical wave functions 

[(V.l)(x)b-q(x) b ] *({*/}) = 0. (2.4) 

where 

[ V 3 l A x )] b = hi x )b - Kj{x - ja) h c lj(x - ja) c . (2.5) 
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3 The axial gauge; splitting the Hamiltonian 

Next we discuss the axial-gauge-fixing procedure. This is most easily done on a cylinder 
of dimensions L 1 x L 2 , with open boundary conditions in the x 1 -direction and periodic 
boundary conditions in the x 2 -direction pQ. For any function /(rc^x 2 ), we require 
that /(x 1 ,^ 2 + L 2 ) = f(x 1 ,x 2 ). The range of coordinates is x 1 = 0, a, 2a, . . . , L 1 , 
x 2 = 0, a, 2a, . . . , L 2 — a. For any physical state Gauss's law implies that 

h(x l ,x 2 )^ = [JZ^x 1 - a,x 2 )h(x l - a,x 2 ) - iV 2 l 2 )(x l - a, x 2 ) + q(x)] * , 

so the operators on each side of this expression may be identified. Taking the gauge 
condition U\(x l ,x 2 ) = 11, which is possible everywhere on a cylinder, we sum over the 
1-coordinate to obtain 



/ 1 (x 1 ,x 2 ) = ^g(y 1 , a ; 2 )-^(P 2 / 2 )(y 1 ,; 

3/!=0 y!=0 



(3.1) 



which is the lattice analogue of (jl.ip . There is a remnant of Gauss's law not determined 
by (|3.1|) . which is the condition 



x 



Lr^o 



x 1 =0 



* = 



(3.2) 



on physical states The condition (J3.2|) is the lattice analogue of (jl.2j) . 

In the axial gauge, using the nonlocal expression (J3.1|) for the electric field in the 
x 1 -direction (which henceforth will be called the horizontal direction), the Hamiltonian 
becomes H = Hq + Hi, where 



Hn 



L 2 -a 
x 2 =0 



o[x 



(3.3) 



with 



x 1 ^ 



L Y -a 



H (x 2 ) = J2 ^Hx\x 2 )} 2 -J2T1- [^U 2 {x)^U 2 {x l + a,x 2 ) + c.c] 

x=0 Z9 a 



(3.4) 



and 



Hi 



Ml! 

2a 



L 2 - 



E E 

x 2 =0 x 1 ^ 1 ^ 



X 



x [l 2 (x , x 2 ) — TZ 2 (x , x 2 — a)l 2 (x , x 2 — a) — q(x , x 2 )~\ 
x [l 2 (y\ x 2 ) - n 2 {y\ x 2 - a)l 2 (y\x 2 - a) - q{y\ x 2 )} 



(3.5) 



and where we have now introduced the second dimensionless coupling constant g' , 
defined by (e') 2 = (g' ) 2 /a. The only interaction in the x 2 -direction (which from now 
on will be called the vertical direction) is due to Hi. 

If g' vanishes, the cylindrical lattice splits apart as shown in Figure 1. The dashed 
line at the top of the unsplit lattice on the left in this figure indicates that this line is 
identified with the bottom line. No such identification is made in the split lattice on 
the right. 



x 



L 2 



2a 







a 2a L 1 



x 



Figure 1. The splitting of the lattice at g' Q — 0. 

The operators H (x 2 ) are Hamiltonians of SU(2) x SU(2) principal chiral nonlinear 
sigma models, as noted in reference pQ (see also reference [3]). Each sigma model is 
represented by a horizontal ladder of plaquettes on the right-hand side of Figure 1. We 
see that there is a ladder for each value of x 2 . Setting g' to zero results in decoupled 
layers of (l+l)-dimensional sigma models. Increasing g f leads to an interaction between 
the vertically-separated layers. This fact was used to give a set of simple arguments for 
confinement for g' Q -C g . The expression for Hi, acting on physical states, with q = 0, 
is identical to the expression on the right-hand side of equation (3.6) in reference pQ, 
by virtue of ()3.2j) . 

For readers who are not familiar with Hamiltonian strong-coupling expansions, we 
remark that they start by neglecting the second term of ()2.2|) or ()3.4|) . This term is 
reintroduced by Rayleigh-Schrodinger perturbation theory, yielding the strong coupling 
expansion in 1/<7q. 
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In a strong-coupling expansion, the magnetic flux is allowed to flow through space 
unconstrained. By including the "magnetic" or "plaquette" term of order 1/g 2 , pertur- 
batively, this unphysical assumption is corrected for. In contrast, we start by allowing 
the 1-component of the electric field to flow through space unconstrained. By including 
H\ (we shall discuss how, at the end of the next section) we correct for this assumption. 

At this point, the reader has noticed that our aim is to exploit a type of dimensional 
reduction from (2 + 1) to (1 + 1) dimensions. This reduction is very different from 
"compactification" , that is, making L 2 small. It is in some sense the opposite of Fu 
and Nielsen's idea of a "layer phase" JH] , i n which the coupling between lattice layers 
is made strong instead of weak. The reduction is similar to the "deconstruction" of 
Arkani-Hamed, Cohen and Georgi JH] i n that the difference between the (1 + 1)- 
dimensional Hamiltonians and the (2 + l)-dimensional Hamiltonians in that the sigma 
models can be regarded as gauged and coupled together through the external gauge 
field. 

4 The effective action for the electrostatic potential 

Before we can make much use of the axial-gauge formulation, we need to examine H 
and Hi in more detail. First we consider H Q . If we adopt the interaction representation, 
time derivatives of operators 21 are given by <9 2l = i[H , 21]. By working out the time 
derivative of U2(x 1 ,x 2 ), we find 

l(x\ x 2 ) b = ^ Tr t b d U(x\x 2 ) U(x\ x 2 )^ , 
% 

K{x\x 2 ) b c l{x 1 ,x 2 ) c = ^ r lTt b U{x\x 2 ) ] d Q U{x\x 2 ) , (4.1) 

9o 

where we have dropped the subscript 2, since there is only one spatial component of 
the gauge field. The time dependence of operators is implicit in these expressions. The 
(1 + l)-dimensional SU(2) x SU(2) principal chiral sigma model of the field U G SU(2) 
has the Lagrangian 

£pcsM = 7 ^?TTr^f/ t ^f/, //,*> = 0,1. (4.2) 

The left-handed and right-handed currents are, respectively, 

fi(x) b = Ym h dJJ{x) U( X y , j*(x) b = iTr t b U( X yd,U(x) . (4.3) 
The Hamiltonian obtained from ()4.2j) is 

#pcsm = J dx^{[ti{x) b f + [j 1 L (x) 6 ] 2 } = j dx^{\jf{x) b f + \jf{x) b ] 2 }. (4.4) 
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By comparing (|4.1j) with ()4.3|) . we can see that .ffpcsM i n ()4.4|) is the continuum limit 
of H (x 2 ) in (jS3J). 

Next we turn to Hi. This interaction Hamiltonian is nonlocal, but can be made 
local by reintroducing the temporal component of the gauge field. In one continuous 
infinite dimension R, the function g(x x — y 1 ) = {x 1 — y l \/2 is the Green's function of 
the "Laplacian" ; in other words 

-d 2 1 g(x 1 -y 1 )=5(x 1 -y 1 ). 

On our lattice, with x 1 and y 1 taking values 0, a, 2a, . . . , L 1 , the same function 
g(x l — y 1 ) = \x l — y 1 \/2 is the Green's function of an (L 1 /a + l)-dimensional operator 
An.a, by which it is meant 

t' ! 

In the continuum limit a — ► and thermodynamic limit L 1 — > oo, Ax,i a — > — c^. We 
use this operator to introduce an auxiliary field ^(a; 1 , a; 2 )^ to replace ()3.5|) by 



i 2 =Oi'=0 ^ 

- (g'o) 2 [hix 1 , x 2 ) - 1Z 2 {x l ,x 2 - a)/ 2 (x 1 , x 2 - a) 

- g(x 1 ,x 2 )]$(x 1 ,x 2 )} . (4.5) 

Let us assume that there are only two color charges - a quark with charge q at site u 
and another quark with charge q' at site v (note: the gauge group is SU(2), so it makes 
no difference if we have a pair of heavy quarks or a heavy quark and antiquark). For 
small lattice spacing, we approximate the sum over x l only as an integral to obtain 

,1 ™2\ 



Hi = J2 I ^ 1 ^^9 1 $(x 1 ,x 2 )9 1 $(x 1 ,a; 

L 2 —a „ 

J2 / ^ i [jo l (^ 1 ^ 2 )*(^ 1 ^ 2 )-joV^ 2 )$(^ 1 ^ 2 +«)] 

„9> r, >/ 



2 L 2 -a 
9 ' ~ 



9 °' X *=0 

+ (g' ) 2 qMu\ u 2 ) b - (^o) V^Cv 1 , v 2 ) b . (4.6) 



We wish to stress that we are not really taking a — > 0, but only assuming that a 
is small. Though the expression (|4.5jl makes the physical meaning of the interaction 
clearer than ()4.5|) . we shall keep the regulator, at least implicitly. From the coupling to 
charges, we see that is proportional to the temporal component of the gauge field 
A Qb . We will call $b the electrostatic potential for this reason. 
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From ()4.6|) we see that the left-handed charge of the sigma model at x 2 is coupled 
to the electrostatic potential at x 2 . The right-handed charge of the sigma model is 
coupled to the electrostatic potential at x 2 + a. 

We now state the mathematical problem we wish to solve. In the presence of a 
quark at u and an antiquark at v, what is the effective action ©($), after integrating 
out U7 If u 2 = v 2 , that is, the quarks are only separated horizontally, the effective 
action is given by 

e i©(*) = (o\Te-^ dx ° Hl |0> , (4.7) 

where the state |0) is the tensor product of sigma- model vacua (in other words, it is 
the vacuum of Hq) and where T denotes time (x°) ordering. If u 2 7^ v 2 , the expression 
()4.7j) is no longer correct. In that case, the expectation value needs to be taken with 
respect to different eigenstates of H , as we shall discuss in Section 7. 

The effective action may be expanded in terms of vacuum expectation values of 
products of currents: 

L 2 — a „ 

i e ($) = -i / d 2 x(g' ) 2 a 2 d 1 ^(x° } x\x 2 )d 1 ^(x°,x\x 2 ) 

- E/ d " x j d 2 y[(0\Tri(x°,x\x 2 ) b ri(y ,y\x 2 ) c \0) 

+ (0| T Jo R (x°, x\ x 2 ) b j*(y°, y\ x\ |0>] $(*°, x\x 2 ) b <S>{y\ y\ x 2 ) c + 0($ 4 ) 

dx° [{g',) 2 q{x°)Mx\u\u 2 ) b - (g' ) 2 q' (x°) b <S>(x° ,v' ,u 2 ) b ] 

+ i5 , wzwN(Q') + i5'wzwN(9') > (4-8) 

where d 2 x = dx°dx 1 , d 2 y = dy°dy 1 and Swzwn^) is the Wess-Zumino-Witten-Novikov 
action of a single SU(2) quark charge (see for example references (201) ^ ne details of 
which are not important for our purposes. We will determine the two-point correlators 
of currents in ()4.8|) at large separations. 

If we ignore the quantum corrections in ()4.8|) . the potential between the quark and 
the antiquark is 

V{u l - v 1 ) =a\u l -v 1 \, = q 2 W ^ - 3 ^ 



2 



(g' ) 2 a 2 2 a 

This is the result of Section 6 of reference pp. In the next section, we show the quan- 
tum corrections from the current- current correlators will drastically change this result. 
Physically, these correlators correspond to transverse (that is, vertical) fluctuations of 
the electric string, as we discuss in Section 7. 
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5 The leading-order corrections to the effective ac- 
tion 



From the exact result for the two-point form factor (jA.35|) . discussed in the appendix, 
we will determine the correlation functions 

D(x, y) bc = (0| T Jo L (x°, x\ x 2 ) b j L (y°, y\ x\ |0> 

= (0| T Jo R (x°, x\ x 2 ) b j*(y°, y\ x% |0> , (5.1) 

in (j4.8j) for the SU(2)xSU(2) ~ 0(4) sigma model. We will then use this result to find 
the string tension for horizontally-separated color charges. 

The exact correlation functions will also have contributions from two- as well as all 
higher-point form factors. The complete formula for the Wightman (non-time-ordered) 
expectation value of two operators in terms of form factors is 

m{x)C{y) |0) = (0|«8(o;)|0>(0|£( 2 /) |0> 
d9\ ■ ■ ■ cLOm 



M=l 



(2vr) M M! 



8(aO \9 m Jm-- -^ujiWxj!, . . .,9 M ,j M \ €(y) |0).(5.2) 



To obtain this result, we used the resolution of the identity (|A.3|) . Time-ordered 
expectation values can be written in terms of Wightman functions. In a field theory 
with a mass gap m, the largest contribution to ()5.2j) for large m\x — y\ comes from 
the terms with the smallest number of particle exchanges M. For our problem, with 
two-sigma model charge densities, the first term, i.e. the vacuum channel, gives no 
contribution; therefore we can consider just the case of M = 2. We will evaluate (J5.1j) 
in this way. Viewing the two-point form factor as a vertex between the electrostatic 
potential $ and two excitations of the sigma model, called Faddeev-Zamolodchikov or 
FZ particles, we consider the one-loop diagram: 




This diagram is infrared and ultraviolet finite. A scale is set by the sigma-model mass 
gap m. We will expand this diagram in derivatives of $. Lowest order in derivatives 
means lowest order in momentum, in the Fourier transform of this amplitude. The 
two-particle form factors should therefore be sufficient. If we wanted many higher- 
derivative terms in the effective action of the electrostatic potential <5($), this would 
no longer be the case. 

Taking the complex conjugate of the expression for the form factor ()A.36j) . in the 
appendix, with x replaced by y, applying (|5.2j) . truncating M > 2, and finally ordering 
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in the time coordinates x° and y°, yields the following for the current-current correlation 
functions (|5.1|) 

D(x,y) bc = ^—^ J de lC ie 2 (cosh^ - cosh# 2 ) 2 \F{0 2 - ^)| 2 

x exp {— imsgn(x° — y°) [(x° — y°)(cosh9i + cosh# 2 ) 

- (x 1 - y 1 ) (sinh^i + sinh# 2 )]} , (5.3) 

where sgn(x° — y°) is defined as 1 if x° > y° and — 1 if x° < y°. 

The integration in (j5.3J) can be made somewhat easier by introducing new param- 
eters Q = (0! + 2 )/2 and u = {Ox - 6 2 )/2: 

D(x ,y)bc = o~~ f dVLdu sinh 2 Vt sinh 2 10 

1X 2 J 

x exp{— 2im coshf2sgn(x° — y°) [(x° — y°) cosher — (x 1 — y 1 ) sinhu;]} 

x exp — / ttf 1 — coshf cos . (5.4) 

io ( COSh 2 I V ?r / 

The action ©(<&) is nonlocal, but the excitations of the sigma model are massive 
particles. Hence we expect that ©(<&) is dominated by local terms, obtained from 
the derivative expansion. Let us introduce the new coordinates X^ and r M by x M = 
X M + r^/2 and y^ = — r M /2. We expand <3>(x) and $(y) in powers of r M , in the 
standard way as 

$(x) = $(X) + y^ t $(X) + ^d^iX) + ■■■ , 

<%) = $(X) - y^$(X) + ^<9 M ^$(X) ± ■ ■ ■ , (5.5) 

where now denotes d/dX^. 

Taking care to sum over L and R, the term we want to evaluate in the effective 
action, which is quadratic in the fields, is given by 



ie«(.) = -(g 4 E/**M* + 5.*-5) h 



x9(X + -,*-)9(X --,*•). (5.6) 

At the risk of overemphasizing a point, we remark that the truncation of ()5.2j) to 
M = 2 and the use of the derivative expansion ()5.5|) have the same justification. Both 
are valid approximations in a massive theory. Though we expand in r in (|5.5|) . it is 
not a short-distance expansion in the usual sense. We can integrate over r precisely 
because large m\r\ contributions are suppressed. The result of substituting ()5.5|) into 
(|5.fj|) is really a small-momentum expansion. 
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After substituting (|5.4jl and ()5.5|) into ()5.6|) we integrate over r M . The integrals over 
r M are are of the form 

I{Q ,Q U A) = J (fre-'^^+^^^Air^r 1 ) , 

for some polynomial A(r°,r r ), where 

Q = 2m cosh Q cosh u , Q± = 2m cosh Q sinh uj . 

We therefore need to evaluate I{Qo, Qi,A), for a few choices of A(r°, r 1 ). 

Suppose that ^(r ^ 1 ) has no dependence on r 1 . Then the r 1 -integration will 
produce a term proportional to 5(Qi), which, in turn, is proportional to 5 (sinh a;). 
From the factor sinh 2 co> in the integral expression (|5.4j) . this will give no contribution 
to in lETfl) . Since we are working to quadratic order in r° and r 1 , we therefore 

need only consider A^^r 1 ) = r 1 ,r°r 1 , (r 1 ) 2 . 

It is elementary to show that I(Qo, Qi, r 1 ) = 0. 

The integral I(Qo, Qi, r°r 1 ) is proportional to 5'(smh.to). Again, the factor sinh 2 u; 
in the integral expression (|5.4jl insures that this will give no contribution to i©( 2 )($) 
in (JSU). 

The only integral remaining is 

l[Qo, Qu (r 1 ) 2 ] = - 2 -p^ = -P.V.^0 + 2,5(Q )5"( Ql ) . 
Only the principal value contributes to the effective action, which is 

i©( 2 )($) = -?™ f^V V / I ^ Sinh ^ Sinh u 2 " ^(2mcosh^sinhc,) 
71 VW ^7 i cosh ft cosh w 

x exp- f ^ 6 fl-cosh£cos— ^) fafcpf, a; 2 ) 6 l 2 

7 £ COsh 2 § V 7T / 

+ higher derivative terms . 
The integration over Q and u; can now be done. After some work, we find 



i©( 2 )($) 



1 / % 



3m 2 7r 2 \go 

L 2 -a 

X 

c 2 =0 



1 



exp 



_ 2 / ^e"«tanh 2 
Jo Z 



/ d2x [dMX,x 2 ) b ] 2 + h. d. t., (5.7) 



which is the central result of this paper. 

The correction to S($) which is cubic in $ can be shown to vanish by symmetry 
considerations, 6^ 3 ^($) = 0. The quartic correction &^(&), does not vanish. This 
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can also be obtained using form-factor methods, though the calculation will be longer 
than that above. An interesting aspect of (5( 4 )($) is that it couples fields together 
at neighboring values of the vertical coordinate, e.g. x 2 and x 2 + a. It may lead to 
interesting nonlinear dynamics of electric strings. At this order in g' the form factors 
and the mass spectrum will be altered [2T] . which will also need to be implemented. 



6 The horizontal string tension 

With our result (|5.7|) serving as the second term of the expression (|4.8jl for the effective 
action we find 

L 2 —a „ 

&(^) = -K^2 / d 2 x [di$(x, x 2 ) b ] 2 + h. d. t. + 0($ 4 ) 

9 rx J 



X 2 =0 



dx° [(g' ) 2 q(x\^x ,u\u 2 ) b -(g' ) 2 q\x\^x\v\u 2 ) b ] 

+ <Swzwn(<7) + SwzwnCqO j (6-1) 
Where the factor K is given by 



-2 / — e~ 5 tanh 



(6.2) 



Notice that there is no induced mass term in $. This is not hard to understand; it is 
due to the fact that the FZ particles of the principal chiral sigma model are adjoint 
charges, hence do not screen quarks. 

Notice that at leading order, there are no time derivatives of $ in ©($). The 
effective Hamiltonian is 

L 2 -a „ 

E = / dxl [di$(x,x 2 ) b ] 2 + h. d. t. + 0($ 4 ) 

+ (gtfq&iu 1 , u 2 ) b - (g'ofq&iv 1 , u 2 ) b . (6.3) 
Thus the horizontal string tension is 



°H — 7W&J 



4^' K 8K ' 

where we used, as before, q 2 = (q') 2 = | for our normalization of SU(2) charges. 

Let us look a bit more closely at the factor K. In the asymptotically-free SU(2) x 
SU(2) nonlinear sigma model, the mass of the FZ particles depends on the coupling g 

as 

m=%^e- 2 ^ + ■••), 
a 
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where C is a non- universal constant, which depends on the cut-off method. Thus the 
second term in K according to (|6.2j) is significant. Our result for the horizontal string 
tension is therefore 



n 3 ( 9 '° 

2 V a 



4 0.7296 (gg 47r/3 g 
1 + 3 C 2 tt 2 gl 



(6.4) 



Notice that for very small go, the exponential dominates the denominator, even if 
g' <C go- Thus ()6.4|) depends on the couplings as 



** - ,5onl • (6-5) 



(0.7296)8 a 

This expression is unusual in that all ^-dependence has disappeared. 

If we rewrite our expression for the string tension ()6.4|) in terms of the continuum 
couplings e = goy/o- and e' = g' y/a, we see an expression which is different than (J1.3jh 
In terms of these constants and the lattice spacing, that 



3 (e 



f\2 



aR= 2 



4 0.7296 (ef 4T/(A) 
3 CH 2 e 2 



' 9CV e 2 _ 47r/(e2a) 



(0.7296)8 a 



where the approximation is valid for small a. We cannot take a continuum limit of our 
string tension, holding e and e' fixed. The reason the naive argument leading to (jl.3j) 
fails is in the assumption of no ultraviolet divergences. What our method shows is that, 
at least for the anisotropic range of couplings we consider, there are such divergences. 
The resulting dependence of the string tension on the couplings is not analytic. 

We believe that ()6.4j) cannot be extended to the isotropic regime g' ~ go. To see 
why, imagine that we generalize the regularized theory to one with three couplings. 
In the Euclidean Wilson lattice gauge theory, with lattice spacing a and link fields 
U^(x) G SU(2), where x lies on a three-dimensional lattice, this means an action of the 
form 

S = J2 jV Tr U^U^x + p.a)Up(x + ua)U v {x) , (6.7) 

where /t denotes the unit vector in the /i-direction and g^ u = g Ufl . There are three 
distinct couplings in this model. The regime analogous to that we consider is got 
902 = 9\i- There are other regimes, we could apply our result ()6.4|) to, namely goi <^ 
<7oi = 9vi and g\i go\ = go2- Clearly there must be a crossover phenomenon between 
different behaviors of the string tension. It seems plausible that there is other crossover 
behavior between these regimes to g i ~ go2 ~ 9i2 in which (|1.3|) takes place. Since 
the string tension does not depend on g' to our order of approximation, it may be that 
the crossover occurs at a value of g' /go which is not extremely small. 
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7 Physical aspects of excitations 



In this section, we complete the picture of the confining phase by discussing the vertical 
string and the nature of the pure-glue excitations. 

We already presented the basic mechanisms of linear potential between vertically- 
separated quarks and the area decay of space-like Wilson loops in reference p. We 
will show how these mechanisms fit into a general picture of the excitations. 

Let us begin by asking what the excitations are if go > and g' = 0. As we 
discussed in Section 3, the system splits up into decoupled layers of nonlinear sigma 
models. The possible excitations are massive FZ particles which can travel horizontally, 
but not vertically. There is, however a restriction on the these excitations, which is 
that the residual gauge-invariance condition (|H.2jl must be satisfied. If we approximate 
the sum over x 1 as an integral, this condition states that for each x 2 

y dx 1 [jfo\a?) b -j*(x\x l -a) b ] - g 2 Q(x 2 ) b ^ = , (7.1) 

where Q(x 2 ) b is the total color charge from quarks at x 2 and is any physical state. If 
there are no quarks, the total right-handed charge of FZ particles in the sigma model 
at x 2 — a is equal to the total left-handed charge of FZ particles in the sigma model at 
x 2 . 

Let us picture two-dimensional space partitioned into a set of parallel horizontal 
layers. Each layer contains a sigma model. The FZ particles move within a layer, as 
in the left-hand side of Figure 2. Now suppose we increase g' from zero to a small 
value. Since an FZ particles at x 2 has left-SU(2) charge at x 2 and right-SU(2) charge 
at x 2 + a, horizontal electric strings must join the FZ particles together. The strings 
lie between the layers. Because the constraint ()7.1|) is satisfied, these strings can be 
consistently introduced. We now have a similar picture of excitations, but now with 
strings with the tension calculated in the last section. This is shown in the right-hand 
side of Figure 2. We see now that the vertical electric flux is carried by the FZ particles 
themselves - they are short segments of vertical electric flux. In fact, if we introduce 
a quark and antiquark with a vertical separation, there will simply be a line of FZ 
particles and horizontal strings joining them together. 

The term in the effective action of the electrostatic potential we calculated in Section 
5 is due to charge fluctuations in the sigma model. According to the picture we have 
just outlined, this means it is due to transverse (that is, vertical) fluctuations of the 
string joining a quark-antiquark pair. 

The reader should not be misled by the right side of Figure 2 into thinking that 
there is no 1-component of the electric field if g' Q = 0. There is electric field produced 
by the FZ particles in this case, but this field carries no energy. 
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Figure 2. The circles represent FZ particles in the layers between the dashed 
lines. On the right, electric strings connect the FZ particles. 

Naively, the vertical string tension is simply pQ 

a v = - « — 2 e- 2 ^ . (7.2) 
a g a z 

This result follows from assuming that the energy in each layer is m. Strictly speaking, 
the right-hand side of (|7.2|) is a lower bound to the vertical string tension, ay > m/a. 
We expect that there will be corrections to this formula. These can be found, in 
principle, just as we did as for the horizontal string tension. There is a subtle difference 
however. In the expression for the effective action. Instead of (|4.7j) . we must calculate 



e 



iS($) = lim (m\Te- l ti dx ° H ' . (7.3) 



The state $ is not the vacuum. We must now have an FZ particle in every layer 
between the quark and antiquark. This problem is under study. 



8 Conclusions 

By spitting (2 + l)-dimensional SU(2) Yang-Mills theory anisotropically into integrable 
models, with an interaction between these models, we obtained an expression for the 
string tension (in one direction only) at small values of the couplings. 

Two problems yet to be solved for this gauge theory are finding a more precise result 
for the vertical string tension (discussed in the last section) and the mass gap of the 
gauge theory. The latter problem is not straightforward, as the glueball excitations are 
collections of many FZ particles. Furthermore, we do not expect the number of these 
FZ particles is fixed. The right approach may be to obtain a better understanding of 
the vacuum state first, then examine gauge-invariant correlators in this state. 



17 



The effective action we used to find the horizontal string tension contains further 
nonlinear terms, as we mentioned at the end of Section 5. The quartic term contains 
mixed correlators of left-handed and right-handed charges of the sigma model. It is 
necessary to use (|A.38|) . as well as (jA.36|) There is no reason why this term cannot be 
determined; all that is needed is sufficient effort. Perhaps this term could shed light on 
the crossover phenomenon discussed at the end of Section 6. It will also be important 
to include corrections to the mass spectrum and the form factors at this order fU\ . 

We could extend our methods to SU(iV) gauge theories, if we knew the form factors 
for the SU(iV)x SU(iV) sigma model with N > 2. Our basic idea can be formulated for 
any value of N pQ. The form-factor problem for N > 2 is tough, in part because of the 
presence of bound states of the fundamental FZ particles. The S-matrix was worked 
out some time ago j22j . and agrees with the Bethe Ansatz approach for an equivalent 
Fermionic model j2H|. The S-matrix becomes unity in the 't Hooft limit iV — > oo, g$N 
fixed, and the form factors should simplify in this limit (note: there is a large- N limit of 
the SU(iV)x SU(iV) sigma model, whose S-matrix is non-trivial [21]. This is a different 
model, because the limit is not the standard 't Hooft limit). 

The picture of confinement and excitations described in Section 7 suggests that the 
non-Abelian gauge theory may be dual to another field theory. The weak coupling 
diagrams would correspond to the strong-coupling terms of the dual theory. 

A splitting similar to the one we have used can be done in (3 + 1) dimensions. There 
is an important difference, however. As in (2+1) dimensions, some electric-field compo- 
nents squared are included in the interaction Hamiltonian. The new feature is that the 
interaction Hamiltonian also contains a magnetic-field component squared. Our meth- 
ods would therefore not yield, strictly speaking, a weak-coupling result. Investigations 
in this direction may be of some value, nonetheless. 
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Appendix: Exact S-matrices and form factors 

We now present an introduction to how S-matrices and form factors are exactly de- 
termined for the O(A^) nonlinear sigma model in (1 + 1) dimensions. Though not a 
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comprehensive treatment of (1 + l)-dimensional S-matrices and form factors, this re- 
view is self-contained. Since most people who might find this paper of interest are not 
experts on these subjects, we felt it was necessary to provide complete explanations of 
their less obvious aspects, which are not abundantly available. We hope this appendix 
is sufficiently complete for the reader to understand the results used in the remainder 
of the paper. 

Integrability of the (1 + l)-dimensional quantized asymptotically-free O(N) non- 
linear classical sigma model became of interest, after it was first established by 
Pohlmeyer for the classical case Then general arguments were made that the 

infinite set of dynamical charges can be generalized to the quantized theory [27J . The 
generation of a mass gap in the 1/iV-expansion was discovered considerably earlier |2E]. 

The Lagrangian of the O(N) sigma model depends on an iV-component vector field 
field n = (n 1 , . . . , n ), of unit length n ■ n = 1: 

C N LS=7T2' r r^n-d v n. (A.l) 

For N = 4, the identification between (jA.lj) and the SU(2) xSU(2) Lagrangian is 
made by U = n 4 ll — in b cr b . 

From the 1/iV-expansion, we know that there are N basic species of massive particle, 
which we label by letters j, k,l,m, . . . taking the values 1, 2, . . . , N. The particle states 
are eigenstates of momentum q,p, . . . as well as species j,k, . . . , created on the vacuum 
by Faddeev-Zamolodchikov operators or FZ operators 21(g) t: 

\q,j,p,k,...)=%(q)]Vl(p)l--\0) . 

We are using the Heisenberg representation, so we work with in-states and in-operators 
or out-states and out-operators. 

From either the the 1/iV-expansion or from the assumption of integrability, we find 
that particles are neither created nor destroyed when scattered (that is, there is no 
particle production). Hence multi-particle scattering may be decomposed as a sequence 
of two-particle scatterings. The requirement that this decomposition is consistent is 
the Yang-Baxter equation or factorization equation. For nonintegrable theories, there 
is no such decomposition. 

We can find expressions for Green's functions which are valid at large distances us- 
ing exact form factors, for models with a mass gap. Green's functions are not exactly 
known, except those for the spins of the Ising model, where all the form factors are 
known [22] • Nonetheless, an expression can be written down for the vacuum expec- 
tation value of a product of operators which approaches the correct answer at large 
separations. 
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Al. The exact S-matrix of the O(iV) nonlinear sigma model 



The two-particle S-matrix should have, on fairly general intuitive grounds, the 
following form: 



where the four-tensor in species indices S(s), depends on the center-of-mass energy 
squared s = (q + p) 2 , and ± refers to Bose or Fermi statistics (we are going to consider 
the former only). We introduce the other Mandelstam variable t = (p' — p) 2 and 
the rapidities 61,62, related to the momenta by the standard relation p = mcosh#i, 
Pi = msinh^i, q = m cosh 62, q\ = msinh# 2 . The relative rapidity is (9 12 = 62 — 61- 

Bound states of two particles can form only when the center-of-mass energy is less 
than the total rest-mass energy. Above the the total rest-mass energy, the energy 
spectrum is continuous. Thus there may be poles in S(s) in the complex s-plane for 
s < Am 2 , and there is a cut in S(s) for s > Am 2 . Assuming crossing s ^ t, there is 
a cut for t > Am 2 . Kinematically, s = 2m 2 (l + cosh# 12 ) and t = 2m 2 (l — cosh# 12 ). 
Hence there are really two cuts in the complex s-plane; one for s > Am 2 and another 
for s < 0. A minimal analyticity assumption is that any poles present lie on the real 
s-axis on the interval < s < Am 2 . 

If we make the change of variable from s to 6 = 6yi, we find the analyticity structure 
is as follows. The s-plane is mapped to the region bounded by Im# = and Imd = ir, 
called the physical strip (the physical strip is not the same as the physical region, which 
is the on-shell kinematically allowed region of 6). Each of the boundaries Im^ = 0, ir, 
—00 < Re# < 00 is a cut. Poles lie on the Re# = axis only. There may be poles 
on this axis outside the physical strip, but these do not correspond to physical bound 
states (there are no bound states for the model we are studying). Note that crossing 
corresponds to 6 — ► in — 6. The resolution of the identity or overcompleteness relation 
in rapidity space, which is straightforwardly obtained from that in momentum space, 
is 



H = |0)<0| + X; / 1 • " " \6 M)]M . . . , 9 l ,j 1 )(e ll j 1 , . . . , 6m , 3m I . (A.3) 



The tensor S(6) can be decomposed in the following way, following Zamolodchikov 
and Zamolodchikov [TH], |3Uj : 



out 



(q',j,p',l\q,m,P, k) in = 8(q[ -pi)5(p[ - qi)S^ k 

±5{q[-q l )5{p' 1 -p 1 )St 




(A.2) 



M=l 




SL{0) = S^SxiO) + S l m 6iS 2 (6) + sid l k s 3 



(6). 



(A.4) 



We represent this decomposition pictorially as 
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s» k (9)=s 1 (e)X + s 2 {6) / x x + s 3 {9)\<^ = 

km km k m k m 

Each line connecting two indices in this picture is a Kronecker delta, in accordance 
with diagrammatic lore. The rapidity 9 is drawn as an angle between incoming and 
outgoing lines (though its range is, of course, not restricted from to n). 

The functions Si(9), S 2 (9) and S 3 (9) are assumed to be real on the real #-axis. 
Therefore, by the Schwartz reflection principle, Si(— 9) = Si(9)*, S 2 (—9) = S 2 (9)* and 
S 3 (-9) = S 3 (9)*. 

If the reader stares at the diagrammatic form for S(9), he or she should be able to 
see that under crossing, S 2 (in — 9) = S 2 (9) and 5*3 (in — 9) = S\(9). 

There are two non-trivial conditions satisfied by the S-matrix elements which, to- 
gether with maximal analyticity, will give its complete determination. The first of 
these is unitarity within the two-particle sector. This is satisfied because there is no 
particle production. Applying unitarity to (jA.2|) yields 



St k {-9)Sf{9) = 5ffi . 

This expression can be worked out by multiplying Kronecker deltas together and con- 
tracting indices, or by diagrammatic methods. In either case, we obtain the three 
relations 

S 2 {9)S 2 {-9) + S 3 {9)S 3 {-9) = 1 , (A.5) 
S 2 (9)S 3 (-9) + S 3 (9)S 2 (-9) = , (A.6) 
NS^S^-O) + S 1 (9)S 2 (-9) + S 1 (9)S 3 (-9)+S 2 (9)S 1 (-9) 

+S 3 {9)S 1 {-9) = 0. (A.7) 

The second nontrival condition is scattering factorization or the Yang-Baxter relation. 
This condition makes the three- and higher-particle scattering amplitudes unique. If 
we imagine three classical particles scattering, there are two possible orderings in which 
this scattering can occur. Each of these orderings corresponds to a certain decompo- 
sition of three-particle S-matrix elements in terms of two-particle S-matrix elements 
in the quantum theory. If either decomposition is correct, they both are; we therefore 
identify them. 

Let us imagine three particles scattering with particle rapidities 9\, 9 2 and 9 3 . The 
relative rapidities are 9 — 9 2 — 9±, 9' — 9 3 — 9 2 and 9 + 9' = 9 3 — 9\. The Yang-Baxter 
equation is shown pictorially in Figure 3. We have labeled each of the incoming particles 
with a species index i, j, k and each of the outgoing particles by another species index 
/, m, o. 
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Now comes the real work - turning the Yang-Baxter equation into something man- 
ageable. Each side of the Yang-Baxter equation can be decomposed into twenty-seven 
terms, each proportional to a six-index tensor, which is a product of three Kronecker 
deltas. The left-hand and right-hand sides of the Yang-Baxter equation are decom- 
posed pictorially in Figures 4 and 5, respectively. Each of the terms is equal to a 
product of three functions S a (0) times a product of three Kronecker deltas. We abbre- 
viate products of three such functions, e.g. S2(9)S^(9 + 6')S2{0') as S2S3S2, keeping 
the arguments in the order 9,9 + 9', 9' always. Though there are twenty-seven terms 
on each side, each term must be one of only fifteen distinct tensors, up to the factor 
of three S's. That is because there are exactly fifteen ways to make a six-index tensor 
from a product of three Kronecker deltas. 




Figure 3. The factorization equation for three-particle S-matrix elements. 

For N = 2, there is an additional subtlety, which is that the fifteen tensors con- 
structed from products of three Kronecker deltas are not linearly independent. This 
is important for understanding the sine-Gordon/massive- Thirring model, but is of no 
relevance to this paper. These tensors are linearly independent for iV > 3. 

With the labeling of species as in Figure 2, we can see the form of each term in 
Figures 4 and 5. For example, the first term on the left-hand side of the Yang-Baxter 
equation is seen from Figure 4 to be Si(9)S\(9 + 9')S'i{9')5ij5 k 5 mo . In this way, the 
Yang-Baxter equation reduces to fifteen algebraic equations. Of these fifteen equations, 
the seven which are proportional to one of the following: 

8\5f8 k , 5 l °5™8 ik , 8 l k 5™5°, 5 l k 5 m °5 l3 , 6 lm 6°6, k , 5^5°, 5 l k 5™5° , 

are trivially satisfied. Of those eight remaining, five are redundant, leaving three non- 
trivial equations [HI], [HOI- The terms proportional to ^5^5° and each give 

S2S3S3 + S3S3S2 = S3S2S3 . (A. 8) 
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^ \J ^ 




Figure 4. Expansion of the left-hand side of the factorization equation. 
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Figure 5. Expansion of the right-hand side of the factorization equation. 
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The terms proportional to 8 lo 5™8j]- and S lm SikS° each give 

S2S1S1 + S3S2S1 = S3S1S2 • (A-9) 

The terms proportional to S l jS mo S lk and 6 lo 6™6ij, also give (|A.9J) . but with the arguments 
9 and 9' reversed. Finally, the terms proportional to 5\5 mo 5jk and S S^S^ each give 

NS1S3S1 + S1S3S2 + S1S3S3 + S1S2S1 + S2S3S1 + S3S3S1 + SiS%Si = S3 Si S3, (A. 10) 

In each of ()A.8|) . (jA.9|) and ()A.10|) the arguments are 6,8 + 6' and respectively. 

Next, let us solve the unitarity and factorization equations. If the function h(6) is 
defined as h(6) = S 2 (6)/S 3 (6), $EM becomes 

h(9) + h{9') = h{6 + 6') . 

Unless h{6) vanishes, the only possible solution is 

h{6)= l -6 ) (A.ll) 

for some constant A. Note that ()A.6|) yields h(—6) + h(6) = 0, which is automatically 
satisfied. Defining another function p{6) = Si(6)/S 3 (6), equation (jA.9|) and (jA.ll|) 
imply that 

\8 

P {6+6r l +p{er l = -j- 

The solution for p{9) is 

m = — ^ , (A.12) 

IK — 9 

where k is another constant. In terms of the original functions, 

Ss(&) = ^S 2 (e) , S 1 (9) = -r^S 2 (9) . (A.13) 

(7 IK — U 

Next we substitute (jA.13|) into ()A.10|) and multiply both sides by 9{ik -9)(9 + 9'){ik - 
9 — 9')9'{ik — 9'). The result is that a fourth-order polynomial in rapidities 9 and 9' is 
zero. The zeroth-, first- and fourth-order terms are identically zero. Both the second- 
and third-order terms are zero provided 

X(N-2) 

By the crossing property and (|A.13j) . we have that k = it and thus A = 2n / (N — 2). 
By (lA~5l) and (TQ31) . we obtain 

S2(9)S 2 (-9)= 9 \ n2 , 

° "T (N-2) 2 

= -W^e^ e) • SM = - { N-l)i-e) S * W ' (A ' 14) 
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Using crossing, we can write the first of these equations as 



s 2 {e)s 2 {m + e) = - i 4?r2 . (A.i5) 

U + (7V-2) 2 

We will solve (IA.15J) for the solution of the two-particle S-matrix inside the physical 
strip, assuming maximal analyticity. Rather than following references ^B] and [301 at 
this stage, we will instead use an elegant prescription invented by Karowski, Thun, 
Truong and Weisz for the sine-Gordon/massive-Thirring model |31j . We use this pre- 
scription not only because it is straightforward, but because it gives the solution in 
precisely the form we need for obtaining form factors of current operators. We are not 
aware of this method being used for models other than sine-Gordon in the literature 
(we suspect it lies buried in the notes of the serious practitioners of the subject), but 
the result is certainly well known, and can be obtained by other methods. 

To use the prescription of reference [3T], it is convenient to decompose S(9) differ- 
ently. Instead of (|A.4j) . we write 

S{6) = P S {6) + P S S S (9) + P A S A (6) , (A.16) 

where Pq, P$ and P A project to the singlet, symmetric-traceless and antisymmetric 
irreducible representations, respectively, on N <g> N : 

(Poy mk = ^ l 5 mk , {Ps) j L = \{5{5 m -5iJ k ) , 
{Pa)L = \ {S 3 k S l m + 8151) - ±5H mk . (A.17) 
The equations (jA.14j) . (jA.15|) become 

s A (e)s A («i + e) = (i + £)(i+ iA N 92 



iA 



e j v e + Tri j e 2 + a 2 1 - 



iA ' 



m- — {e + i x ){ 9-m) — SA{dh Ss{d) ~ Wi\ SA{dh (A - 18) 

where, as before, A = 2ir/(N — 2). We will use the first of (|A.18|) to solve for S A (9) . 
Following reference [31] we assume 

1. The function S A (9) is analytic 1 and non-zero in the interior of the physical strip 
(it has, of course, cuts on the boundaries). 



1 In the case of the sine-Gordon model, there are bound-state poles on the imaginary axis in the 
analogous S-matrix element, for some choices of the coupling (22, |30| . |34j . Hence this analyticity 
assumption is only valid for the case where the solitons have repulsive interaction. After the answer 
is found for this case, it can be generally applied by analytic continuation 
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2. | In S A (9)/ sinh(z — 9)\ —>■ as |Re z\ — > oo, for any fixed choice of 9 in the physical 
strip. 

Now sinh(z — 9) has only one zero of z in the physical strip, provided 9 is in the 
physical strip. Let C be the counter-clockwise contour enclosing the physical strip; so 
C extends from — oo to oo and from 7ri + oo to 7ri — oo. Then 

1hSa(0) = ^ / . J Z . lnS A (z) 
2tti J c smh{z — 9) 



1 f°° dz 



2iri J _ oa sinh(z 



\n[S A (z)S A (z + 7ri)] 



1 /' x dz 1+ iA 



2vri sinh(z - 0) 1 - £ 



In— (A-19) 



where in the last step, we used (IA.18J) . 

Though ()A.19|) is a succinct expression for the S-matrix, it is not yet the most useful 
form. To find a better form, we need the following Fourier transforms 

dz p i?z ipi(0-*-V2) r°° dz f \\\ 1 

' ^e^lnfl--) =-(\-e^ x ) , 



, 2tt sinh(,2 -9) 2 cosh f ' J-00 ^ V 2 / t 

which can be worked out using basic complex-integration methods. Substituting these 
Fourier transforms into (jA.19|) . we obtain 

S A {9) = exp j 6 '^' 1 sinh (g) . (A.20) 

This is the form of the S-matrix we shall use to discuss form factors. We note that 
this expression (jA.20j) can be converted to the Zamolodchikovs' rational expression of 
gamma functions JH], pO] using the integral formula 



T(z) = exp / -J- 



e -tz 



+ (z-l)e 



Rez > 



which can be checked by differentiation and comparison with the integral formula for 
the psi function and verification that the right-hand side is equal to one at z — 1. This 
integral formula was first used by Weisz e 36, to write the infinite-gamma-function- 
product-ratio expressions of Zamolodchikov [32], [HD! for the sine-Gordon S-matrix in 
a more compact form. 

We have obtained the minimal 0(A/")-symmetric S-matrix, i.e. that with as much 
analyticity as possible. When expanded in powers oil/N, this agrees with the S-matrix 
of the nonlinear sigma model obtained by standard 1/A^-expansion methods |37j - 
If the two-particle S-matrix elements are multiplied by CDD factors [HBj 

t-i- sinh 9 + i sin 
sinh 9 — i sin ctk ' 

k 



27 



unitarity and factorization are unaffected. One such S-matrix obtained in this way 
is that of the O(N) Gross-Neveu model. See references ^H], [HI] for more discussion. 
The supersymmetric sigma model is a nonlinear sigma model and a Gross-Neveu model 
coupled together; its S-matrix was obtained in reference |39|. 

A2. Exact form factors 

The determination of exact two-particle form factors was initiated by Vergeles and 
Gryanik j^U] for the fundamental particles of the sinh-Gordon model and by Weisz 
plo] for the solitons of the sine-Gordon model. Subsequently, Karowski and Wiesz |3j 
obtained generalizations of this result for the sine-Gordon model and extended the 
method to other models. Smirnov and later Kirillov and Smirnov 41^ found extensions 
to higher-point form factors. Smirnov formulated a set of axioms underlying the entire 
subject 02] and which was useful in studying specific field theories. In the meantime, 
Smirnov's axioms have actually been proved as theorems, assuming maximal analyticity 
and the LSZ axioms jlHj. Though establishing their validity from deeper principles is 
certainly worthwhile, the validity of Smirnov's axioms can be argued from symmetries, 
crossing, integrability and a little physical intuition. 

Our interest in form factors is that they contain off-shell information about inte- 
grable quantum field theories. This means that they can be used to study deformations 
of such theories which are no longer integrable. The form-factor program has lead to 
nonperturbative calculations in statistical mechanics [29j, 0^ an d condensed-matter 
physics jlH], which agree well with experimental measurements. 

A form factor is a matrix element of an operator 05 (x) between multi-particle states. 
We can obtain all of these from 

. . . , e M ) h ... JM = (o| Q5(o) \e M ,j M , 0i,ji)- m (A.21) 

by crossing and the Lorentz-transformation properties of the operator 05 (x). For ex- 
ample 

M 

. . . , 6 M ) n ... jM exp -iJ^Pii ■ x i = (°l \0mJm, • • • , di, Ji} in (A.22) 
1=1 

Let us state some basic commutation relations of the FZ particle creation in- 
operators and annihilation in-operators ^(0)^ and 2l J (6 l ), respectively (the arguments 
of these operators are now rapidities instead of momenta): 

%(02)%(0i) =Sl?(6 2 -6 1 )% n (6 1 )%(6 2 ) , 

%(02)%(0iY = 2tt%5(0 2 - 00 + Sf(9 2 - l )K m (6 i y%(6 2 ) . (A.23) 
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The third equation is obeyed for a free field theory with the S-matrix element replaced 
by the identity on two-particle space. The 27r in the normalization of the first term on 
the right-hand side of this equation follows from the overcompleteness relation (jA.3|) 
in rapidity space. The form of the second term on the right-hand side of this equation 
follows from crossing. 

For pedagogical reasons, we shall list the form-factor axioms for general multi- 
particle form factors and attempt to convince the reader that they are reasonable. Our 
list follows Essler and Konik [IS], though we attempt to provide more justification for 
some of the axioms. We use only four of these five axioms. We will only consider the 
case of operators with mutually-local commutation relations. For discussion of how to 
deal with other commutation relations, see Chapter 6 of Smirnov's book [12] or the 
article by Essler and Konik [46J. 

1. Scattering Axiom. This follows from the properties (jA.23|) of the FZ operators. It is 
/ ^(#1, • • • , Oi-i, Oi+i, Oi, 6i+2 ■ ■ ■ , 9M)h-k-ik+i3ik+i"3M 

= ^jik+i ~~ ^+l)/ >B (^l ) • • • ; @i-U 0t> 0i+l 3 #i+2, • • • , (>m )j : -j, :'.',/.•..,./...,-./.•,;• (A.24) 

This axiom is sometimes called Watson's theorem [TFj . 

2. Periodicity Axiom. This axiom is subtle. It is motivated by crossing. The axiom 
states 

f*(0i, ■ ■ ■ , ( 'm)j:-j-,, = f*{0M ~ 27ri, 6>i . . . , QM-i)j M j v ~jM-i ■ (A.25) 

To see where this axiom comes from, let us see what happens when a creation FZ 
operator in front of a ket (state vector) is replaced by an annihilation operator behind 
a bra (dual state vector) by crossing. Consider the Green's function of FZ operators 
and 23 

(o| %(0O ®(o) % M {e M )% M .A0M-^ ■ • •%(^ 2 ) t |o) connected 
= <o| %M 23(0) % m (0m)% m ^m-iV ■ ■ -%MY |o> 

- (O|%(0 1 )®(O)|O)(O|% f (^M) t %, f _ 1 (eM-i) t ---%(^) t |O) , 

which is "connected" in the sense that the vacuum intermediate channel is subtracted 
[23] ■ This Green's function can be thought of as M — 1 incoming particles being 
absorbed by a vertex corresponding to the operator 03(0) which then emits a single 
particle. Consider the pair of particles, with labels 1 (the outgoing particle) and M. 
Under crossing, these both become incoming particles, but with 9\ replaced by 9\ + 7ri. 
To see this, notice that this change preserves all the relativistic invariants (pj ±pj + i) 2 , 
j = 2, . . . , M — 1, but it switches the two invariants S\m = (pi + Pm) 2 and t 1M = 
(Pi ~Pm) 2 - Thus 



(0\%M) ®(0) % m (0m)% m ^(0m-J ■ ■ ■ |0) connected 

> @M)ji—j M 



/*(#! + vri A,..., M.Wm ■ (A-26) 
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Suppose that instead of switching the invariants Sxm and £xm> we switch the invariants 
S12 = (pi + P2) 2 and t X2 = (pi - p 2 ) 2 - Then we find 

(0| 03(0) 2l iM (^) t %- M - 1 (^-i) t • • •% 2 (^ 2 ) t |0) c 



f*(9 2 ,9, Om-0. ■ ~i) ;: ,.. nw: . (A.27) 



' connected 

Comparing (|A.26|) and (|A.27|i . the axiom follows. 

3. Annihilation Pole Axiom. This axiom is important for relating form factors of M 
particles to form factors of M — 2 particles. Though we will not apply M > 2 form 
factors in this paper, we mention the axiom anyway, since higher-point form factors 
should be eventually be of interest in our problem. The idea is the following: let us 
again consider the form factor f^(9i, . . . , 0M)ji-j M - Now there is the possibility of the 
(M — l) st and M th particles annihilating (they can be antiparticles) at 9m = @m-i + 71 "! 
(so that sm-im — {pm-i + Pm) 2 will vanish). Thus there must be a pole in the form 
factor which corresponds to this annihilation at 8m = Om-i + 7ri. There are two terms 
in the residue of this pole: i) The (M — l) st particle may scatter with the (M — 2) nd , 
(M — 3) rd , . . . , 1 st particles before annihilating the M th particle, or ii) the (M — l) st 
particle may not scatter with any particles before annihilating the M th particle. The 
axiom is 

i Res/ 35 (6'i, . . . , 0M)ji-j M \ dM=eM _ 1+n i = /^(^l) •••■> ®M-2)h"3M-i C 3M-i3M 
- b hj 1 iPl - VM-l)\ h {82 - Vm-1) ■ ■ ■ b tM _ 3jM _ 3 [0M-3 ~ Om-1) 
X S jM-ijM-2^ AI - 2 ~ ^M-l)/ 3S (^l) • • • > &M-2)ki-k M -2 Ck M -ij M > (A.28) 

where C is the charge-conjugation matrix. The normalization of the left-hand side 
follows from the standard state normalization, e.g. < 6\6' >= 2tt5(6 — 9'). 

4. Lorentz-Invariance Axiom. If an operator 05 (x) carries Lorentz spin s, the form 
factors must transform under a boost 9j — > 9j + a for all j — 1, . . . , M as 

f™(9i + a,...,9 M + a) = e sa / !8 (0i, ...,9 M ), (A.29) 

We hope that the reader will have no trouble distinguishing spin s from center-of-mass 
energy squared s. 

5. Minimality Axiom. Just as we assume S-matrix elements have as much analyticity as 
possible, so we make a similar assumption of form factors. In order to check the validity 
of this principle, all that can be done is to compare with some perturbative method, 
which means either standard covariant perturbation theory or the 1/A-expansion. In 
cases where we can find form factors, minimality stands up very well. If a first guess 
for the form factor f^(9i, . . . , 9m) satisfies the first four axioms, then so does 

f® (n n \ - f®<7) n ^({cosh^- - 9 k )}) 

/minimal^, ■ ■ ■ Jit) ~ f (*!,..., ( {cOS h(^. - }) ' (A - 30) 
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where Pm and Qm are symmetric polynomials. In this way, we can eliminate all the 
poles in form factors, except those corresponding to bound states. In order for Axiom 
3 to be satisfied: 



PM\8 M= M _ 1+ ri - Pm-2 , QM\e M= e M _ 1+7T i - Qm-2 ■ 

The overall normalization of the form factors is not determined by these axioms 
The normalization can be found for the case of current operators, which is what we 
shall apply to the (2 + 1) -dimensional Yang-Mills theory. 

A3. Currents of the O(N) nonlinear sigma model 

After this review of form-factor concepts, we next apply these ideas to the case of 
two-particle form factors, which we shall find explicitly for currents of the OA) sigma 
model. The two-particle form factor may be written as 



/ (0ij02)jua = e s 1 F (6 2 - 0i) jlj2 , 

by Axiom 4. In terms of the function F^(9i, O^)^-, Axioms 1 and 2 are 

F*{e) klk2 =e-* e S{^ 2 {6)F*{-6) nn , 
F*(2m-6) nn =e se F*(6) nn , (A.31) 

respectively. 

Let us now recall the projectors Pq, Ps and Pa defined in (|A.17|) . We have already 
used the fact that these diagonalize the S-matrix acting on species tensors transforming 
as N <g> N. We define 

which allows us to rewrite (|A.31|) as 

FfoiO) = e s9 F*(2ni - 6) , F*{0) = -e s *Ff (2th - 9) . (A.32) 

We can solve the equations f)A.32|) for the form factors, up to an overall normalization 
by a contour- integration method 

Consider a function F{9) satisfying 

F(27ri - 6) = -e s9 F{6) , F{9) = e sf> S A {6)F(-6) . 

Define the contour C to be that from — oo to oo and from oo + 27ri to — oo + 27ri. Then 

i r^m f dz , z-9 , , f°° dz , z-9 , F(z) 
In F{6) = — : coth In F(z) = I — coth In ■ 



c 



Am 2 w 7_ oc 4vri 2 F(z + 2ni) 
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Differentiating this formula with respect to 9 (we will explain why in a moment) 



dF{9) f 30 dz 1 , F(z) 

In 



d9 8vri sinh 2 ^ -e s8 F(-z) 

00 dz 1 
oo 8tti smh 2 ^ 

From our expression for SU in (jA.20|) and differentiating the integral formula 

00 dz , z-0 , sin 2 %i-#) 1 

— coth sinh — = ^— — — (A.33) 



47ri 2 ni sinh£ 2 sinh £ 



(which can be done using basic complex-integration methods) with respect to 9, and 
finally integrating dF(9)/d9 with respect to 9, 



-2£ 



d£ e^- 2 — 1 sin 2 ^ 



7T1 



F(e)=Gexp2 X f^n — £hf— • (A - 34) 

where G is a constant. The reason we had to differentiate and integrate with respect 
to 9 is that otherwise, the second term of (|A.33|) will not lead to a convergent answer. 

Now we consider the form factors of the O(N) current operator of the sigma model. 
This operator is J^ fc = n- ? 9 M n fc — n k d fl n J . By Hermiticity, translation invariance, anti- 
symmetry and Lorentz invariance 



<0| J**(0) \9 2 ,m,9 h l)=iG(8l8 k m - <&tf)(p 



1 -P2, 



f°° dE e— - 1 sin 2 f (ni - 9 12 ) , , . 

xexp2 l t^ti — ikr^- (A - 35) 

This is the expression for the current-operator form factor with as much analyticity as 
possible. The normalization of the right-hand side of (jA.35|) is obtained by the crossing 
relation (jA.26|) for the two-point case. We will fix this normalization for the case of 
N — 4, which is the case we wish to study further. It is not difficult, however, to obtain 
the normalization for any N using slightly different methods. 

The normalization is fixed by examining the matrix elements of a charge operator. 
The eigenvalues of this operator are fixed by symmetry considerations. In this way the 
value of the form factor can be specified for particular rapidities. We will find that 
G = l. 

As we already mentioned, the connection between the SU(2)-valued field U and the 
unit four- vector n is U = n 4 — inVj,. The relation between the currents defined in ()4.3|) 
for the principal chiral sigma model and those for the vector sigma model is therefore 
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The left-handed charge density obeys the algebra 

bo (*% jo(y%] = iv^e^^x 1 - y^jiix 1 ), , 

and the left-handed charge is Q b = J dx 1 Jo L (x 1 ) b . The charge therefore obeys the 
commutation relations 

[Ql Q L C ] = iV2e bcd Q L d , 

Since the charge is in an orbital representation, the possible eigenvalues of Q\ are 
(isospin zero), 0, ±\/2 (isospin one), etc. We find the same relations for the right- 
handed charge. 

From (jA~22j) and (jA~35j) . we have 

(0| jo' R (x) b \62J2, 0i,ji) = iV2G (Sj^Sjzb - S hA S jlb ± e bjlj2 )m(cosh9i - cosh0 2 ) 
xexp{-im[x o (cosh0i + cosh0 2 ) - x 1 (sinh 9t + sink 9 2 )]} F(9 2 - 9 X ) , (A.36) 

where the plus or minus sign corresponds to the left-handed (L) or right-handed (R) 
current, respectively, and 

pvm 9 r d£ e-* - 1 sin 2 fdf r« . 2 e(7ri-g) /A 

F# =exp2 / — 7 — = exp — / — sin (A. 37 

KJ io ( ^ + 1 sinh£ P J £ cosh 2 | 2vr 1 ; 

Now under crossing, the relation (jA.26|) for our simple two-point case yields 

<^i,ii|j'o' (^)&|^2,i2> = iG\/2 (5 jl4 5 hb - S h J jlb ± e bjlh ) m(cosh6 l i + cosh6> 2 ) 
x exp {— im[x° (cosh 9\ — cosh# 2 ) — x 1 (sinh^ 1 — sinh# 2 )]} 
xF{9 2 -m + 9 1 ) , (A.38) 

Integrating each side of (jA.38|) over x 1 yields 

(0i » Ji| Qb' R 1^2, J2> = 1GV2 (S hi S j2b - 8 hi 5 hb ± e bjlj2 ) (9x1 6 2 ) ■ 

This result yields the eigenvalues of Q b ' R for isospin one, provided G = 1. 

Equation (|A.35|) with G = 1 is a result of Karowski and Weisz who also checked 
its validity in the 1/A-expansion. We shall use it to obtain an expression for the 
effective action of the electrostatic potential $. 

When an interaction is added to the action of an integrable model, the form factors 
will recieve corrections |2Ij. The mass spectrum will be altered as well. We shall 
not work to high enough order in g' Q for these effects to be taken into account. If 
higher-order calculations can be done, however, they will need to be included. 



33 



References 

P. Orland, Phys. Rev. D71 (2005) 054503. 
[2] S. Mandelstam, Bull. Am. Phys. Soc. 22 (1977) 541; Phys. Rev. D19 (1979) 2391. 
[3] M. Karowski and P. Weisz, Nucl. Phys. B139 (1978) 455. 

[4] P.A. Griffin, in Fourth International Workshop on Light-Cone Quanti- 
zation and Non-Perturbative Dynamics, Polana Zgorzelisko, Poland 1994. 
Published in Polana Light-Cone Quantiz. (1994) 240, hep-th/9410243. 

[5] R.P Feynman, Nucl. Phys. B188 (1981) 479. 

[6] P. Orland and G.W. Semenoff, Nucl. Phys. B576 (2000) 627. 

[7] P. Orland, lhep^th/9607134] (1996). 

D. Karabali and V.P Nair, Nucl. Phys. B464 (1996) 135; Phys. Lett. B379 
(1996) 141; D. Karabali, C. Kim and V.P. Nair B524 (1998) 661; Phys. Lett. 
B434 (1998) 103; Nucl. Phys. B566 (2000) 331, Phys. Rev. D64 (2001) 025011. 

[9] P. Orland, Phys. Rev. D70 (2004) 045014. 

[10] J.P Greensite, Nucl. Phys. B166 (1980) 113. 



[11] R.G. Leigh, D. Minic and A. Yelnikov, hep-th/0604060 (2006). 
[12] M.J. Teper, Phys. Rev. D59 (1998) 014512. 



[13] M. Caselle, P. Grinza and N. Magnoli, hep-th/0607014 (2006). 

[14] J. Greensite, Prog. Part. Nucl. Phys. 51 (2003) 1. 

[15] P. Orland, hep-th/0608067 (2006), submitted to Physical Review D. 

[16] A.B. Zamolodchikov and Al. B. Zamolodchikov, Nucl. Phys. B133 (1978) 525. 

[17] M. Creutz, Quarks, Gluons and Lattices, Cambridge University Press (1983), 
Chapter 15. 

[18] Y.K. Fu and H.B. Nielsen, Nucl. Phys. B236 (1984) 167; Nucl. Phys. B254 (1985) 
127. 

[19] N. Arkani-Hamed, A.G. Cohen and H. Georgi, Phys. Rev. Lett. 86 (2001) 4557. 



34 



[20] A. P. Balachandran, G. Marmo, B.S. Skagerstam and A. Stern, Gauge Symme- 
tries and Fibre Bundles: Applications to Particle Dynamics, Lect. Notes 
Phys. 188 (1983), 1-140, Springer- Verlag; A.M. Polyakov, Mod. Phys. Lett. A3 
(1988) 325; H.B. Nielsen and D. Rohrlich, Nucl. Phys. B299 (1988) 471; T. Ja- 
cobson, Phys. Lett. B216 (1989) 150; P. Orland, Int. J. Mod. Phys. A4 (1989) 
3615. 

G. Delfino, G. Mussardo and P. Simonetti, Nucl. Phys. B473(1996) 469. 

E. Abdalla, M.C.B. Abdalla and A. Lima-Santos, Phys. Lett. B140 (1984) 71; 
PB. Wiegmann; Phys. Lett. B142 (1984) 173. 

A.M. Polyakov and PB. Wiegmann, Phys. Lett. B131 (1983) 121; PB. Wieg- 
mann, Phys. Lett. B141 (1984) 217. 

V.A. Fateev, V.A. Kazakov and P.B. Wiegmann, Nucl. Phys. B424 (1994) 505. 

A. M. Polyakov, Phys. Lett. B59 (1975) 87; E. Brezin, J. Zinn- Justin and J.C. Lc 
Guillou, Phys, Rev. D14 (1976) 2615; E. Brezin, J. Zinn- Justin, Phys Rev. B14 
(1976) 3110; W.A. Bardeen, B.W. Lee and R.E. Schrock, Phys. Rev. D14 (1976) 
985; Y.Y. Goldschmidt and E. Witten, Phys. Lett. 91B (1980) 392. 

K. Pohlmeyer, Commun. Math. Phys. 46 (1976) 207; M. Liischer and K. 
Pohlmeyer, Nucl. Phys. B137 (1978) 46. 

I.Ya Arafeva, P.P. Kulish, E.R. Nisimov and S.J. Pacheva, preprint LOMI-E-1- 
1978 (1977); A.M. Polyakov, Phys. Lett. B72 (1977) 224. 

H. E. Stanley, Phys. Rev. 176 (1968) 718. 

B. M. McCoy, C.A. Tracy and T.T. Wu, Phys. Rev. Lett. 38 (1977) 793; M. Sato, 
T. Miwa and M. Jimbo, report Kyoto RIMS-225(1977); B. Berg, M. Karowski and 
P. Weisz, Phys. Rev. D19 (1979) 2477. 

A.B. Zamolodchikov and Al.B. Zamolodchikov, Ann. Phys. 120 (1979) 253. 
M. Karowski, H.J. Thun, T.T. Truong and P. Weisz, Phys. Lett. B67 (1977) 321. 

A. B. Zamolodchikov, Commun. Math. Phys. 55 (1977) 183. 
M. Karowski and H.J. Thun, Nucl. Phys. B130 (1977) 295. 

B. Berg, in International Colloquium of Complex Analysis, Microlocal 
Calculus and Relativistic Quantum Field Theory, Sept. 3-13, 1979 Les 
Houches, France (1980). 



35 



[35] E.T. Whittaker and G.N. Watson, A Course of Modern Analysis, Cambridge 
University Press, Chapter 12, pg. 249. 

P. Weisz, Phys. Lett. B67 (1977) 179. 

B. Berg, M. Karowski, V. Kurak and P. Weisz, Phys. Lett. B76 (1978) 502. 

L. Castillejo, R.H. Dalitz and F.J. Dyson, Phys. Rev. 101 (1956) 453. 

R. Shankar and E. Witten, Phys. Rev. D17 (1978) 2134. 

S. Vergeles and V. Gryanik, Sov. Journ. Nucl. Phys. 23 (1976) 704. 

F.A. Smirnov, preprint LOMI-E-3-1986 (1986); J. Phys. A19 (1986) L575; A.N. 
Kirillov and F. A. Smirnov, Phys. Lett. B198 (1987) 506; F.A. Smirnov, Int. J. 
Mod. Phys. A9 (1994) 5121. 

F.A. Smirnov, Form Factors in Completely Integrable Models of Quantum 
Field Theory, Adv. Series in Math. Phys. 14, World Scientific (1992). 

H. Babujian, A. Fring, M. Karowski and A. Zapletal, Nucl. Phys. B538 (1999) 
535; H. Babujian and M. Karowski, Nucl. Phys. B620 (2002) 407. 

P. Fonseca and A.B. Zamolodchikov, J. Stat. Phys, 110 (2003) 527. 

M.J. Bhaseen and A.M. Tsvelik, in From Fields to Strings; Circumnavigat- 
ing Theoretical Physics, Ian Kogan memorial volumes, Vol. 1 (2004), pg. 661, 
M. Shifman, A. Vainshtein and J. Wheater ed., cond-mat/0409602 

F.H.L. Essler and R. Konik, ibid., pg. 684, |cond-ma t /0412421 

K.M. Watson, Phys. Rev., 95 (1954) 228. 



36 



